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Topological phase transitions can occur in the dissipative dynamics of a quantum system when 
the ratio of matrix elements for competing transport channels is varied. Here we establish a relation 
between such behavior in a class of non-Hermitian quantum walk problems [M. S. Rudner and L. S. 
Levitov, Phys. Rev. Lett. 102, 065703 (2009)] and nuclear spin pumping in double quantum dots, 
which is mediated by the decay of a spin-blockaded electron triplet state in the presence of spin-orbit 
and hyperfine interactions. The transition occurs when the strength of spin-orbit coupling exceeds 
the strength of the net hyperfine coupling, and results in the complete suppression of nuclear spin 
pumping. Below the transition point, nuclear pumping is accompanied by a strong reduction in 
current due to the presence of non-decaying "dark states" in this regime. Due to its topological 
character, the transition is expected to be robust against dephasing of the electronic degrees of 
freedom. 
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Since the first observation of spin blockade in verti- 
cal GaAs double quantum dots^, spin- blockaded trans- 
port has been observed in a variety of systems such as 
lateral double quantum dots in GaAs, Si, and Si/SiGe 
heterostructures^^— , InAs nanowires^, and carbon nan- 
otubes-. Much of this work was driven by the need to 
better understand the coupled dynamics of electron and 
nuclear spins in double dot systems with potential ap- 
plications in spintronics and quantum computation. A 
variety of interesting and surprising phenomena such as 
current bistabilities and hysteresis^'^'^, very long time 
scale switching^, and periodic oscillations^ have been ob- 
served, and their origins linked to the dynamical polariza- 
tion of nuclear spins (DNP). Although the involvement of 
nuclear spins in these phenomena is clear, in many cases 
the underlying mechanisms remain a mystery. 

Spin blockade of dc transport occurs in a two-electron 
double quantum dot when the electron spins form a 
triplet state that prohibits both electrons from occupy- 
ing the same site. In this case, as shown in FigH^, resid- 
ual current arises from mechanisms that do not conserve 
the electron spin such as the spin-orbit interaction, and 
the hyperfine coupling to nuclear spins in the host lat- 
tice. Because the hyperfine contact interaction conserves 
the total spin of all electrons plus nuclei, each hyperfine- 
mediated electron spin flip is accompanied by a nuclear 
spin flip in the opposite direction. In the presence of 
strong spin-orbit coupling, e.g. as in InAs systems^;-, the 
physics of DNP can be very different from that studied 
previously in the absence of spin-orbit coupling^ii^. In 
particular, by making several transitions between singlet 
and triplet states using a combination of hyperfine and 
spin-orbit processes, the decay of a single electron spin 
can lead to a change of nuclear polarization by an amount 
which can have either sign, and a magnitude potentially 
even greater than one unit of angular momentuiiiii. 

In this paper we explore the rich quantum dynamics 
of coupled electron and nuclear spins that results from 
the coherent competition of hyperfine and spin-orbit de- 
cay channels in spin-blockaded quantum dots. We focus 




FIG. 1: Competition between hyperfine and spin-orbit decay 
in spin-blockaded double quantum dots, a) A triplet state 
decays via a hyperfine or spin-orbit mediated transition to a 
singlet state. The singlet state is coupled to the drain, and 
decays with rate 7. b) Energy levels and transitions. We focus 
on the subspace indicated by the dashed oval, c) Inhomoge- 
neous hyperfine coupling due to non-uniform electron density. 
Red dotted line shows approximation to smooth density pro- 
file consisting of uniform density shells in which nuclear spins 
collectively couple to an electron as separately conserved "gi- 
ant spins." 

on the polarization transferred to the nuclear spin bath 
by the decay of a single electron in one of the blockaded 
triplet states. This quantity reveals a strikingly high sen- 
sitivity of the DNP production efficiency to the presence 
of the competing spin-orbital decay channel. 

In order to study DNP in this regime, we develop a 
class of models which capture the essential physics of po- 
larization transfer during electron spin decay. One of 
the main difficulties in describing nuclear pumping in 
the presence of spin-orbit coupling stems from the lack 
of a concrete conservation law which directly relates the 
changes of electron and nuclear spin polarizations. Any 
model of this process must account for the possibility 
of multiple electron spin transitions which can lead to a 



variety of final nuclear spin states. For a typical system 
containing N sa 10^ nuclear spins, the exponentially large 
Hilbert space makes exact analytical or numerical solu- 
tions difficult to obtain. However, as we show below, the 
problem can be made tractable by introducing approx- 
imations which greatly reduce the number of variables 
while retaining the key degrees of freedom responsible 
for the mechanism of polarization transfer. 

We begin by employing the "giant-spin" model, as used 
e.g. in Ref . |12| . in which the electrons interact with a sin- 
gle large collective spin formed from all the spins in the 
nuclear spin bath. This model describes the case where 
the local hyperfine coupling to each nuclear spin within 
each dot takes on a uniform value, A. Within this approx- 
imation, the problem near the singlet-triplet resonance, 
circled in Fig[TJ), can be viewed as a one-dimensional 
hopping problem in the space of polarization of the gi- 
ant collective nuclear spin, see Figl2^. We study this 
model numerically, and obtain additional insight from 
comparison to the solution of a related quantum walk 
modeU^ in which the quantity analogous to DNP is de- 
scribed by a topological invariant which takes on integer 
values. We find a non-analytic dependence of the polar- 
ization transfer on the ratio of hyperfine and spin-orbit 
coupling strengths, with complete suppression of DNP 
when the spin-orbit coupling exceeds the net transverse 
hyperfine field. This behavior is a direct manifestation 
of the topological phase transition which occurs in the 
quantum walk modeU^. 

To investigate the role of inhomogeneous hyperfine 
coupling, we then employ a model in which the nonuni- 
form hyperfine couplings are approximated by d shells 
of constant coupling, as shown in FigjlJ;. Here, nuclear 
spins couple to form d large collective spins which inter- 
act with the electrons. The resulting dynamics can be 
viewed as a hopping problem in a d-dimensional space 
indexed by the polarizations of each of the d collective 
spins, see FigsHJj and c. 

Although a numerical approach is not possible for 
the general case, exact analytic results for a related 
d-dimensional hopping problem show universal features 
which are independent of the specific grouping of nuclear 
spins. The behavior obtained for rf > 1 is essentially 
analogous to that found in the d — I case (the giant 
spin model). While the details of the behavior near the 
transition are sensitive to the particular decomposition 
into collective spins, the strong suppression of DNP in 
the spin-orbit dominated phase is found to be generic, 
suggesting that it will persist in more realistic models. 
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FIG. 2: Mapping of coupled dynamics of electron and nu- 
clear spins onto quantum walk models describing single par- 
ticle hopping on suitable lattices, a) In the uniform coupling 
model, the state of the conserved "giant-spin" is labeled by 
its z-component m. Spin-orbit and hyperfine transitions with 
amplitudes u and Vm can be interpreted as hopping in a tight- 
binding model defined on a one-dimensional bipartite lattice, 
b) For two groups of nuclei, the nuclear state is labeled by 
the 2-components m\ and m^ of two collective spins I\ and 
/2, and the dynamics are described by a quantum walk on 
a two-dimensional bipartite lattice, c) For an electron den- 
sity profile approximated by d constant density shells, the 
dynamics can be viewed in terms of a quantum walk on a 
d-dimensional bipartite lattice with hopping amplitudes «'"' 
describing hyperfine transitions with each of the a. = 1 . . .d 
collective spins. 

energy when two electrons occupy the right dot. Ener- 
getically, electrons which form a spin singlet can thus as- 
sume either the (1, 1) or (0, 2) charge configuration, with 
one electron on each dot or both electrons on the right 
dot: 
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When the electrons form a spin triplet, the Pauli exclu- 
sion principle forbids double occupation of the right dot. 
Due to tunnel coupling, the singlet states |(1, l)s) and 
1(0, 2)s), Eq.([T]), hybridize to form "bonding" and "anti- 
bonding" states IS") and |5"). In the presence of an ex- 
ternal magnetic field, the triplet splits into its three Zee- 
man sublevels. For concreteness, we consider decay of 
the state \T^) when its energy is close to that of the sin- 
glet states \S), with all other states far away in energy 
(see FigUlD). The relevant electronic states for our model 
are thus 
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I. CONSTRUCTION OF THE MODEL 

We begin by reviewing the relevant two-electron states 
of a spin-blockaded double quantum dot, see e.g. Ref.[l|. 
For simplicity, suppose that the left and right dots each 
support a single orbital state \L) or \R). An applied 
potential bias approximately compensates the charging 



Because of spin-orbit coupling, interdot tunneling is 
accompanied by a spin-rotation which couples the triplet 
and singlet states |T+) and 15*) with an amplitude u. 
In the basis of Eq.®, the purely electronic part of the 
Hamiltonian is written as 



Hn = 
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(3) 



where Et^ and eg are the energies of the states |T_|_) and 
\S), respectively. The imaginary term —ij/2 in is ac- 
counts for the decay of the singlet state due to coupling 
of 1(0, 2)s) to the drain, see FigU^. Without loss of gen- 
erality, we take the spin-orbit coupling matrix element 
u to be real. A microscopic derivation of the value of u 
is beyond the scope of this work. However, because the 
spin rotation occurs during tunneling, u is proportional 
to the admixture of |(0, 2)s) in \S), i.e. to the parameter 
Cq2 in Eq.(I2). In addition, the value of u is sensitive to 
the orientation of the dots relative to the crystallographic 
axes, and to the direction of the applied magnetic field. 
The hyperfine interaction between electron and nuclear 
spins also couples the triplet and singlet states. In a two- 
electron system, the hyperfine Hamiltonian 

i^HF = ^^In • [Sl<5(xi - X„) + S2(5(X2 - X„)] (4) 
n 

couples the spin Si (2) of each electron to each nuclear 
spin I„, with weight proportional to the probability to 
find the electron at the location X„ of nucleus n. Here, 
for simplicity, we consider a single species of nuclear spin, 
but a generalization to multiple species is straightfor- 
ward. The spin- flip terms S^,^)^^ couple the electron 
states with spin projections along the z-axis differing by 
one unit of angular momentum. 

In the basis of Eqs.([2]) and ([2]), Hamiltonian @ takes 
the form 
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(6) 

is the hyperfine coupling weighted by the local electron 
density p(X„) = {ip\S{x — X„)|V') (see FigHJ;), and, due 
to antisymmetry of the wavefunction, ri„ = +1(^1) if 
nucleus n is located in the left (right) dot. Because the 
hyperfine interaction is local, the off-diagonal matrix ele- 
ments of Hamiltonian ^ between |r+) and \S) are pro- 
portional to the amplitude of |(1, 1)5) in \S), i.e. to the 
parameter Cn in Eq.©. 

The sign factors rjn indicate that the difference between 
transverse nuclear polarizations in the left and right dots 
couples the electron triplet and singlet levels. The math- 
ematical annoyance of alternating signs can be removed 
by applying a 7r-rotation about the z-axis to all spins in 
the right dot via the operator U — e^*'^^^^", where the 
sum is taken over all spins in the right dot. In the ro- 
tated frame, the Hamiltonian i?Hp = U'^HyifU takes the 
simpler form: 
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The factor l/\/2 arises from the normalization of the 
states in Eqs.([T]) and ©. Up to these numerical prefac- 
tors, the transformed Hamiltonian is equivalent to that 



describing the hyperfine interaction for a single electron 
in a quantum dot with an electron density profile consis- 
tent with the distribution of couplings {^„}. 



II. GIANT SPIN MODEL (d = 1) 

In the special case where An = A for all n, which cor- 
responds to an electron density that is uniform within 
the dots and zero outside, the square of the total nuclear 
spin operator, I^ — {J2n I")^! commutes with the Hamil- 
tonian. In this case, all nuclei in the system act together 
coherently as one "giant" spin. For fixed /, the configu- 
ration space of the system is then defined by the electron 
states |T+) and |S'), and by the z-projection m of total 
nuclear spin, I^\m) = m|m), with —I < m < I. Com- 
bining Eqs.Q and Q, the Hamiltonian for this system 
can be written as 
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where Ae is the triplet-singlet detuning, and /+(^) is the 
raising (lowering) operator for the giant spin. 

Below we neglect the polarization-dependent Over- 
hauser shift AP by absorbing its mean value into the 
definition of Ae. For polarizations which are not too 
large, this approximation is justified by the fact that the 
typical off-diagonal matrix elements of Hm are of the or- 
der A\/N, while the Overhauser shift only changes by an 
amount of order A when the nuclear polarization changes 
by one unit of angular momentum. For a typical dot con- 
taining N Ki 10^ nuclear spins, the variation of the Over- 
hauser shift thus imposes only a small perturbation on 
the dynamics of the system. In a similar spirit, we also 
ignore the nuclear Zeeman energy, which is assumed to be 
small compared with the inverse lifetime of the blockaded 
state. 

Decay of the blocked triplet state occurs through elec- 
tron spin- flip transitions to the state \S), which is broad- 
ened due to its coupling to the drain lead. These transi- 
tions can be mediated by cither the hyperflne interaction 
or the spin-orbital interaction. The hyperfine process is 
accompanied by a change of the z-projection of nuclear 
spin. Am — ±1, whereas for the spin-orbital process 
Am — 0. As illustrated in Figl2^, the resulting coher- 
ent dynamics in the combined Hilbert space of electron 
and nuclear degrees of freedom can thus be viewed as a 
hopping problem on a one-dimensional bipartite lattice. 

In the basis {\m) ig) \T^/S)}, the state of the system 
IV') is described by the amplitudes V'm = ('7^T'+|V') and 
i/'m = {"fTT- Slip) , and evolves according to the equations of 
motion (with ft = 1) 
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with 



Vm = A^Jl{I + l) -m(m+l). 



(9) 



(10) 



The hopping amphtudes Vm, which originate from the 
transverse hyperfine field, attain a maximum value 
v^i^x = ^\/I{I + 1) for unpolarized states to « 0, and 
become small near maximum polarization \m\ < /, see 
FigHi. 

Suppose the system is initially in the blockaded elec- 
tron spin state, with nuclear polarization tuq. What is 
the average change in nuclear polarization Am = m — mo 
caused by the decay of the electron spin? To formulate 
the problem more precisely, we consider the situation 
where an electron is injected into the triplet state |r+) at 
time i = 0, with an initial nuclear spin state character- 
ized by total angular momentum I and z-projection toq. 
The system then executes a "quantum walk" under the 
equations of motion (jQ]) , with initial state 
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0. 



(11) 



The wave packet describing the quantum walker will 
spread throughout the lattice and leak out through its 
components on the S-sites, decaying completely as i — >■ 
oo. The value of to at the site from which the system 
decays determines the final value of nuclear polarization 
left behind when the electron escapes. Given the proba- 
bility Pjn for the system to decay from each singlet site 
TO., we would like to evaluate 



(Ato) = ^(to - Too) Pra, Prr 



7|V^Wpdt.(12) 



This expression for P™ results from the fact that the non- 
Hermitian equations of motion (J16p lead to decay which is 
a sum over local terms describing decay from each site of 
the lattice, Jj(V'IV') — ^ X^m 7lV'mP- Because the system 
decays completely as i — > oo, ^^ P,„ — 1. 

To explore the behavior of this model, we have solved 
Eq.(j9l) numerically with initial condition (fTTj) for giant 
spins with / = 25, / = 50, and / — 100. The polar- 
ization transfer (Ato,), Eq. (fT2l) . is plotted for each initial 
polarization rriQ in the upper panel of Fig|5)D. In addi- 
tion, we also show the inverse of the average dwell time 
T = — Jq°° t ■^{tp\^) dt for each case in the bottom panel 
of Figlii. 

Two very different situations arise depending on the 
relationship between the spin-orbit coupling matrix el- 
ement u and the maximum hyperfine coupling matrix 
element Umax- In case I, indicated by the dotted line 
u = ui > Umax in FiglHK, spin-orbit coupling dominates 
the dynamics for any initial polarization toq, and nuclear 
spin pumping is strongly suppressed (see Figl^JD). For 
u — uii < Umax, howcvcr, the dynamics can be domi- 
nated either by hyperfine coupling or by spin-orbit cou- 
pling, depending on the value of the initial polarization. 
The system possesses critical points to = ±?7i*, within 



TO* 



±I\ 



(13) 



where, locally, the strengths of hyperfine and spin-orbit 
coupling are nearly equal. For initial polarizations satis- 
fying I Too I < TO,, the dynamics are hyperfine-dominated 




-1 -m, -1/2 1/2 TO* 1 

Initial Condition nio/I 

FIG. 3: Results for the giant spin model, Eq.(|9]). a) Depen- 
dence of hyperfine matrix element Vm on the z-projection of 
total nuclear spin for a giant spin of length I. For moderate 
values of the spin-orbit coupling, u — uu < Umax, the system 
possesses critical points ±r7i*, see Ea. (|13|) . where the spin- 
orbit and hyperfine couplings are roughly equal, b) Average 
polarization change and inverse lifetime versus initial polar- 
ization. For all traces we take Ae — and 7 — Umax. When 
spin-orbit coupling is strong, u = ui = 1.2umax, nuclear spin 
pumping is suppressed. For moderate values of the spin-orbit 
coupling, u = Uu = 2/3umax, (Am) ~ 1 for initial states with 
Umo > uii, and (Am) « for Vmo < uii- Near the critical 
points m — m,/I, the system possesses "dark states" with ex- 
tremely long lifetimes. We include a decay term —171/2 with 
7t ~ 10^''7 in the triplet energy to cut off the divergence of 
the lifetime near the critical points. 

and approximately one unit of angular momentum is 
transferred to the nuclear spin subsystem per electron. 
Outside the critical points, i.e. for |to-o| > to-,, spin-orbit 
coupling dominates and polarization transfer is strongly 
suppressed (see Fig[3lD). As the length of the giant spin 
increases, the distinction between the behaviors in these 
two regimes becomes more sharply defined. In particu- 
lar, (Ato) becomes sharply quantized to 1 for \mo\ < m<,. 
Near the critical points, the dwell time of an electron 
in the system can become very long (see lower panel of 
FigOb). 

The behavior in all of these regimes can be understood 



in terms of a simplified model. If the quantum walk ^ 
only explores a window of sites which is small compared 
with the scale of the giant spin, 2/ + 1, then we may 
approximate the ?7i-dependent hopping amplitudes Vm 
by a single amplitude v = v^g which characterizes the 
strength of the transverse hyperfine field when the giant 
spin has z-projection rriQ. In this same spirit, we also 
extend the lattice to infinity, — oo < m < oo. These 
approximations make the quantum walk translationally 
invariant, and allow us to find an exact analytical solu- 
tion to the dynamics. 

The translationally-invariant model described above is 
identical to the one dimensional non-Hermitian quantum 
walk which was studied in Ref.[13|. In that work, we 
found that the expected displacement (Am) achieved be- 
fore decay is quantized as an integer: 



(Am) 



1 V > u 
u > V . 



(14) 



The quantized value of {Am) is determined by the wind- 
ing of the phase between two components of the Bloch 
eigenstates of the bipartite one-dimensional system as 
the momentum k is taken through the Brillouin zone. 
Equivalently, the value of {Am) can be determined di- 
rectly from the winding number of the complex amplitude 
Ak = u + ve^'' . In the regime where this winding number 
is zero, corresponding to the situation where spin-orbit 
coupling u is stronger than the hyperfine coupling v, the 
result {Am) = indicates that no angular momentum 
is pumped into the nuclear spin subsystem. As we will 
see below, this behavior is quite general and persists for 
more refined models which go beyond the giant spin ap- 
proximation. 

The behavior of the giant spin model, displayed 
in FiglSb, closely resembles the prediction of the 
translationally-invariant model. Some distortions are ob- 
served in the highly polarized regions, |r7i| « /, where 
Vm varies strongly with m,. The striking suppression 
of decay near the upper critical point can be traced 
to the divergence of the lifetime which accompanies the 
topological transition between winding and non-winding 
phases in the translationally invariant model (see RefflSl) . 
The topological transition is manifested in the non- 
translationally invariant system through the presence of 
a topologically-protected "dark" edge state which is lo- 
calized at the phase boundary between winding and non- 
winding phases. This state has zero overlap with the elec- 
tron singlet state, and thus does not decay. Physically, 
the extended lifetime results from the fact that, near the 
critical point, the effective hyperfine and spin-orbit fields 
responsible for electron spin transitions can cancel each 
other. Although suppression of decay is seen near both 
critical points, the effect is much more dramatic near 
the upper critical point where the dark state is stable; 
near the lower critical point, an analogous exponentially- 
growing (delocalized) dark state can be found. 

In addition to the behavior of the lifetime, the small 
"overshoots" in (Am) can also be understood with in- 
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FIG. 4: Phase transition between pumping and non-pumping 
regimes for nuclear polarization described by the giant spin 
model. Average polarization change Am (circles) and lifetime 
f (diamonds) averaged over all initial states — / < mo < / vs. 
the ratio of the hyperfine and spin-orbit al cou plings. The 
black line shows the average displacement (Am), obtained by 
applying the translationally-invariant quantum walk model 
around each initial condition, Eq. (|15|) . 



tuition gained from the translationally-invariant model. 
Semiclassically, the site-dependent hopping amplitudes 
Vm give the "walker" a position-dependent effective mass 
which goes through a minimum at the critical point. As 
a result, the walker experiences a force which attracts it 
to the critical point. For mg < m^, where each electron 
already has a tendency to transfer one unit of angular 
momentum to the nuclear spin subsystem, this additional 
force results in a transfer of more than one unit of an- 
gular momentum per electron, {Am) > 1.^^ Similarly, 
above the critical point, the attraction leads to a nega- 
tive polarization, {Am) < 0. The shape of (Am,) around 
the lower critical point can be understood analogously. 

In experiments, the length / of the giant spin and its 
^-projection mo in the initial state are in general random 
variables picked from the thermal distribution; the length 
/ is distributed according to p{I) ex (2/ -I- 1)^6"^^/^"^ , 
with Iq a constant, while mg is uniformly distributed for 
each choice of /, p{mo) = const. This formula for p(I) 
is obtained for spin-1/2 nuclei using Eq.(l) of Ref.|16| 
and Stirling's formula. In Fig|3]we show the numerically- 
obtained expected polarization transfer (Am) and inverse 
lifetime l/f , averaged over all initial conditions mg, for a 
giant spin with / = 50. The existence of dark states near 
the critical points is manifested in the suppressed cur- 
rent for u < Wmax- Compared with the translationally- 
invariant case, where {Am) displays a sharp step as a 
function of u/(u + v), here the step is rounded into the 
phase where v < fmax- When spin-orbit coupling dom- 
inates, however, i.e. for u > fmax, the suppression of 
nuclear spin pumping is nearly complete. 

The shape of the rounded step can also be understood 
simply within the context of the translationally-invariant 
quantum walk model. For each initial condition mo, we 
take (Am,) = 1 if |mo| < tti* (hyperfine dominates), or 
(Am) = if I mo I > rn^ (spin-orbit dominates). After 



averaging over all initial polarizations ttiq, we find 



m. 



(Am) = Ox (^l--^j+lx-^ = -^, (15) 

where to, is given by Eq. p^ . Expression p5|) is plotted 
as a dotted line in FiglSJj. The good agreement with 
exact numerics for the giant spin model further indi- 
cates that the intuition gained from the translationally- 
invariant model provides a very useful tool for under- 
standing the behavior in the more realistic situation. 



III. BEYOND GIANT SPIN (d > 1) 

Faced with the surprising prediction of complete sup- 
pression of DNP in the spin-orbit dominated phase, it 
is natural to wonder to what extent this result relies on 
the assumption of uniform hyperfine coupling (i.e. on the 
validity of the giant spin model). The giant spin approxi- 
mation tightly constrains the dynamics by truncating the 
dimension of the Hilbert space from 2^ down to 0{Vn), 
where N is the number of nuclear spins in the system. 
This appears to be a rather severe approximation. To ad- 
dress this concern, we now explore a more general class 
of models, derived in a similar spirit, which allow us to 
investigate the effects of nonuniform coupling. 

As illustrated in FiglTJ;, hyperfine coupling is generally 
strong near the center of the dots, where electron den- 
sity is maximal, and weak near the edges, where elec- 
tron density is small. To improve upon the uniform 
coupling model, consider dividing the nuclei into two 
groups, one of "strongly-coupled" spins, and the other of 
"weakly-coupled" spins. In this approximation, depicted 
in FigHb, the nuclear spins within each group form two 
separately conserved collective spins of lengths /i and I2 ■ 
Here the coupled dynamics of electron and nuclear spins 
can be viewed as a hopping model on a two dimensional 
lattice, where the two dimensions index the z-projections 
of the two collective spins, toi and m2. This model cap- 
tures both the transfer of polarization from the electron 
spins to the nuclear spins, and the RKKY-like electron- 
mediated transfer of polarization between the two groups 
of nuclear spins^"— . 

Continuing this reasoning, we can further refine the 
model by approximating the smooth electron density 
profile by d shells of constant density pa, with a = 
1 . . .d (see red dotted line in FigH];, and e.g. Ref.j22|). 
In this case, the nuclear spins couple to form d col- 
lective spins Ii...Id. The corresponding polarization 
transfer dynamics can be viewed as a hopping prob- 
lem on a d-dimensional lattice indexed by the polariza- 
tions TOi , . . . , rud of the d collective spins (see FigllJ;) . 
This "quantum walk" is described by equations of mo- 
tion analogous to Eq.(j9l), with the position variable m 
replaced by a vector m — (to-i, . . . , md): 















where Vm = A-a\/ Ia{Ia + 1) ~ ma{ma + 1) is the hy- 
perfine coupling to collective spin Jq,, with A^ — Ap^, 
and Bq is the unit vector along the axis describing the 
polarization of /q, . 

Our goal is now to calculate the polarization (Am) 
transferred into the nuclear spin bath. 



/•oo 

(Am) = ^(m-mo)P„, P^ = / j\ibi{t)\^ dt,{l 

Jo 
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under the dynamics of Eq. (|16p with initial condition 

i^L^Sm,m„, V4 = 0. (18) 

In particular, we will be interested in determining which 
features of the results are independent of the particular 
grouping into collective spins, and which survive as the 
level of refinement, d, is increased. 

Based on the success of the translationally-invariant 
approximation to the quantum walk in the giant spin 
model (i.e. the case d = 1), we begin by replacing 
the m-dependent hopping amplitudes Vm by constants 
yia) — y^^^ Note that the approximations associated 
with making hopping translationally invariant and with 
extending the lattice of states for each collective spin to 
infinity become more severe as the size of each spin de- 
creases. Thus, although we will proceed for arbitrary di- 
mension d, this number should be considered small com- 
pared to the total number of nuclear spins in the double 
dot. 

The next step is to pass to the momentum represen- 
tation, ^^ = ,^\j § d'^k e^^'^Tp^, where the integral 
is taken over the Brillouin zone — tt < /cq < tt, with 
ae{l,2,...,d}. These Fourier states correspond to co- 
herent nuclear spin states with the transverse component 
of each collective spin a pointing along the azimuthal 
angle k^ . Due to the translational invariance of the sys- 
tem, the equations of motion in the Fourier representa- 
tion break up into 2x2 blocks, one for each momentum 
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dt 






K es 






(19) 



with Ak = u + X)q=i ^ e " . The two-component wave 
functions for different values of k evolve independently, 
and the probability density pk{t) = \ip^{t)\'^ + |V'k(*)P 
to find the system with momentum k at time t decays as 
dtPk — — 7l^k(*)P- '^^^ k-dependence of |Ak| indicates 
that for some giant spin configurations {k^ ~ 0) the ef- 
fective hyperfine and spin-orbit fields add constructively, 
while for other configurations (/cq, ~ tt) they interfere 
destructively. 

Writing TOq, as a derivative with respect to ka via 
■ma ipm = --(^ § '^'^^ di: {e'^'"') V'k and integrating by 
parts to move the derivative onto i/j^, we bring the ex- 
pression for the a-th component of (Am), Eq. p7|) . to the 




1/4 1/2 

u/{u + u' + v") 



FIG. 5: Expected polarization from d = 2 quantum walk 
model, a) Graphical analysis of Eq. (|21[) . For each fixed k\, 
j4k = u + «'e'*i + i)"e'*^ sweeps out a circle as fc2 is varied 
from — TT to TT, encircling the origin for a range around fci « — vr 
bounded by the magenta circles, b) The mean displacement 
(Am2) is equal to the fraction of the Brillouin zone in which, 
for fixed fci, Ak wraps the origin as — tt < ^2 < tt. A sim- 
ilar construction can be used to obtain (Ami) (not shown), 
c) Solution to Eq. (|21|) for (Ami) and {Am2) versus the ra- 
tio of spin-orbit and hyperfine coupling strengths. The ratio 
between the two giant spin hyperfine couplings is fixed to 
the value v' /v" — 2. When spin-orbit coupling it is weak, 
the more strongly coupled giant spin absorbs all spin-flips 
((Ami) = 1, (Am2) = 0). When spin-orbit coupling exceeds 
the maximum hyperfine coupling, u > v' + v" , nuclear spin 
pumping is completely suppressed. 
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(20) 
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The outer integral is taken over the d— 1 momenta kp^a- 
The expression inside the braces is identical to Eq.(5) 
of Ref.[ij| for the displacement in the one-dimensional 
model. As shown there, the value of this integral is 
quantized as either or 1 depending on the winding of 
6k = argjAk} as ka is taken around the Brillouin zone. 

For the one-dimensional case, quantization means that 
the expected change in polarization per electron through 
the system is either 1 if 9k wraps the origin (v > u, 
hyperfine coupling exceeds spin-orbit coupling) or if 
it does not {u > v, spin-orbit coupling exceeds hyperfine 
coupling) . Roughly speaking, the winding of 0k therefore 
distinguishes whether or not the hyperfine coupling is 
strong enough for the electron to flip the nuclear spin. 

To understand the meaning of Eg. (1211) in the multi- 
dimensional case, it is helpful to view the integral in the 
following way: for each fixed configuration of d — 1 col- 
lective spins described by the d — 1 angles {fc^^Q,}, the 
expression inside the braces is either or 1 depending on, 
for the given field of the other spins and strength of spin- 
orbit coupling, whether the electron's hyperfine coupling 
v'"' to the remaining spin is strong enough to induce a 
spin flip (i.e. whether or not 9^ winds the origin as ka 
is varied from — tt to tt, see FigjS^). The integral over 



the remaining d — 1 variables simply counts the "phase 
space" over which this condition is satisfied. This result 
is represented graphically for the case d = 2 in FiglSlD. 



A. Special Case: d = 2 

Recently, the relative dynamics of nuclear spins 
in the two dots of spin-blockaded double quantum 
dots has attracted considerable experimenta P^'^^ and 
theoretica P^i^^" — attention. From a practical point of 
view, understanding the behavior of the difference of po- 
larization between the two dots is important because a) it 
is responsible for dephasing of singlet-triplet qubits, and 
b) because the polarization difference, if carefully con- 
trolled, can be used as a resource to coherently control 
electron spin dynamics^!. The case of our model with 
d = 2, where the electron spins couple to two indepen- 
dent collective spins, is thus particularly interesting if 
the two spins are viewed as representing the nuclear spin 
states in the left and right dots. We analyze this case in 
detail in this subsection. 

The expected displacements (Ami) and (Ar7i2) in the 
two-dimensional quantum walk, see Eq. (PT|) . represent 
the expected amounts of polarization transferred to each 
of the two groups of nuclear spins during the decay of the 
electron spin. The values of (Ami) and (Am2) depend on 
the strength of spin orbit coupling, u, and the strengths of 
the transverse hyperfine fields produced by the two giant 
spins, v' and v" , respectively. Using simple geometric 
arguments based on the construction shown in Fig|5^, 
we find (assuming w' > v") 



(Ami) 



1 

1 - ^ , (Ama) 








v' > u + v" 
\u^v'\ <v" 



(22) 



with 
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,,"2 
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COS 01 



2uv" 



COS0? 



2uv' 



(23) 

Expression (1^^ is plotted in Fig IS}; as a function of 
spin-orbit coupling strength u for the fixed ratio v' /v" — 
2. For strong spin-orbit coupling u > v' + v", nei- 
ther collective spin is pumped at all: we find com- 
plete suppression of DNP in the spin-orbit dominated 
phase just as in the single giant spin model. Interest- 
ingly, for weak spin-orbit coupling we find a quantization 
(Ami) = 1, (A77i2) = 0, which indicates that on average 
the collective spin with stronger hyperfine coupling to the 
electron absorbs the full angular momentum, while the 
more weakly coupled spin is not pumped at all. For an 
asymmetric system composed of one large dot and one 
small dot, electron density and thus average hyperfine 
coupling is larger on the smaller dot. Thus in the limit 
of weak spin-orbit coupling, u <C v' ,v" , nuclear pumping 
may be highly asymmetric, with DNP initially produced 
primarily on the smaller of the two dots. 
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B. General Results 

In contrast to the one-dimensional translationally- 
invariant model, FigjSJ; shows that in higher dimensions 
(Am„) is not strictly quantized as an integer. The 
breakdown of quantization results from the appearance 
of "mixed" phases where the winding number along one 
dimension of the Brillouin zone is 1 for some values of 
the remaining momenta, and for the others. However, 
for the quantum walk (J16l) in any dimension d there is 
always a "non- winding" phase with u > X^a ^ ' where 
all winding numbers are for all values of k in the Bril- 
louin zone. In this phase, (Ama) = for all a. In 
fact, either through graphical methods or with a few lines 
of algebra, one can see that, as u is increased for fixed 
{v^°'^, all of the (Ama) vanish simultaneously at the 
point u = X]q ''^ ■ Thus very generally, strong spin- 
orbit coupling can lead to a dramatic suppression of nu- 
clear spin pumping. 

Furthermore, the expected polarizations (AttIq,) are 
determined purely by geometrical constraints imposed by 
the set of hyperfine matrix elements v^°'' . Importantly, 
just as in the Id case studied in Ref . |13| . both the de- 
tuning e!p — £5 and the decay rate 7 completely drop out 
of the solution. As shown there, the result holds even if 
these quantities are made time-dependent. This implies 
that the suppression of pumping results from nuclear spin 
coherence, and is robust against noise and dephasing of 
the electron singlet and triplet states. 



IV. CONCLUSIONS 

In this work we have analyzed the coupled dynamics 
of electron and nuclear spins in spin blockaded quantum 



dots by mapping it onto a non-Hermitian quantum walk. 
This quantum walk problem can be solved exactly by 
making the approximation that the matrix elements of 
the hyperfine coupling are roughly independent of the 
nuclear polarization m for \m — ttt-qI ^ I, where / is the 
size of a large collective spin formed from a group of nu- 
clei with similar values of the hyperfine coupling to the 
electron. From this solution we find that nuclear spin 
pumping is strongly suppressed whenever the strength 
of spin-orbit coupling exceeds the magnitude of the hy- 
perfine coupling. Numerical simulations show that this 
behavior extends directly to the more realistic situation 
where the polarization-dependence of the hyperfine ma- 
trix elements is included. 



The transition in the nuclear spin pumping efficiency 
is accompanied by an abrupt change in average current 
through the double dot. Due to the presence of "dark 
states" when the hyperfine and spin-orbit processes ex- 
hibit complete destructive interference, the expected life- 
time of the spin-blockaded triplet state is significantly 
longer in the regime where hyperfine and spin-orbit cou- 
plings compete than in a regime where only one of the 
two mechanisms is present. This interference is mediated 
by coherence in the nuclear spin bath, and is therefore 
robust against dephasing of the electronic state. 
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